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Abstract 



J^ ■ We calculate the nuclear relaxation rate l/7i of oxygen in the undoped quasi 

two-dimensional quantum Heisenberg antiferromagnet Sr2Cu02Cl2 above the 
'^ . Neel temperature. The calculation is performed at two-loop order with the 

5 i help of the Dyson-Maleev formulation of the spin-wave expansion, taking all 

Q I scattering processes involving two and three magnons into account. At low 

temperatures T we find l/T^ = ciT^ + C2T^ + 0{T^), and give explicit ex- 
CN . pressions for the coefficients ci (two-magnon scattering) and C2 (three magnon 

f— ^ . scattering). We compare our result with a recent experiment by Thurber et 

^O i al. and show that this experiment directly probes the existence of short- 



wavelength magnons in a two-dimensional antiferromagnet. 
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I. INTRODUCTION 

In a recent experiment [Q Thurber et al. reported nuclear magnetic resonance (NMR) 
measurements of the ^^O and ^'^Cu relaxation rates of the undoped and lightly doped quasi 
two-dimensional quantum Heisenberg antiferromagnet Sr2Cu02Cl2. In this compound the 
anisotropies and inter-layer couplings are particularly weak, so that in the undoped case it 
is an excellent approximation to model the dynamics of the localized S = 1/2 spins at the 
Cu sites above the Neel temperature T/v by an isotropic two-dimensional nearest neighbor 
quantum Heisenberg antiferromagnet. The Hamiltonian is given by 

^='^EE Sr-Sr+a. , (l) 

I" i=x,y 

where the r-sum is over the N sites of the square lattice formed by the Cu-sites, and slx and 
a^ are the primitive lattice vectors with length |aj| = a. The operators Sr are SU{2) spin 
operators satisfying S^ = S[S -\- 1), and J > is the antiferromagnetic exchange coupling. 
NMR experiments measure the relaxation of the nuclear spins due to their coupling to the 
reservoir of the electronic spin lattice. The theoretical framework for the calculation of the 
NMR rates in Heisenberg antiferromagnets has been discussed 30 years ago by Beemann and 
Pincus P], and more recently for two-dimensional antiferromagnets in Refs. 0-|^. However, 
so far only the leading order of the low-temperature behavior has been calculated, which 
is determined by two-magnon scattering and can be obtained from a simple one-loop cal- 
culation. For a quantitative comparison with the experimental result ^ for the NMR rate 
of ^^O it is important to know also the contribution from three-magnon scattering, which 
describes the leading effect of magnon-magnon interactions. In this work we shall explicitly 
calculate this correction and compare it with the experiment of Thurber et al. [Q . 

To begin with, let us briefly derive the NMR rate from Fermi's golden rule. Note that Mila 
and Rice |Q gave the relevant form factor only up to a constant of proportionality. However, 
for a quantitative comparison with experiments we need to know the precise numerical value 
of the form factor. In the simplest approximation, the coupling between the nucleus and the 
electronic spins can be described by an isotropic hyperfine interaction |0-§]. For a nuclear 
spin I^" of ^^Cu at position R this is of the form H'^^ = AcuIr^'Sr, with some characteristic 
energy scale Aqu that can be obtained from experiments. As discussed by Chakravarty and 
Orbach ||^, the ^^Cu NMR-rate is dominated by critical antiferromagnetic fluctuations, and 
cannot be calculated via the perturbative spin-wave expansion. On the other hand, the 
nuclear spins I^^ at the ^^0-site are coupled to the electronic spins of two neighboring Cu 
sites (see Fig. Q), so that the hyperflne interaction is 

^o = Aol^ ■ E SR+b, , (2) 

i=l,2 

where bi = —aix/2 and b2 = aa;/2. Let us emphasize that Eq.(|^) involves the average 
SR+a^/2 + SR_a^/2 of two neighboring electronic spins. Given the fact that the spins couple 
antiferromagnetically, it is clear that at least in the classical picture the total spin vanishes. 
Hence, we expect that the ^^O relaxation rate is dominated by short wavelength antifer- 
romagnetic quantum fluctuations. In this work we shall put this expectation on a solid 



quantitative basis, and show that the NMR relaxation rate l/T^ at the oxygen site probes 
the unique short-distance behavior of magnons in two-dimensional Heisenberg magnets. 

Keeping in mind that 1/Ti is defined such that calculations should always use the matrix 
elements for nuclear spins with magnitude I = 1/2, straightforward application of Fermi's 
golden rule yields for the NMR relaxation rate of ^^O 

1 A^ 

^1 ■^"' af3 i j 

where Ea and |a) are exact eigen-energies and eigen-states of the Heisenberg Hamiltonian 
(0), and Pa = e~^°'^'^/Z, where Z = J2a^~^°'^'^ ■ Here T is the temperature, measured in 
units of energy. In experiments the resonance energy Hun is usually much smaller than 
all characteristic energies of the spin system, so that we may set un = 0. Defining the 
Fourier transformed spin operators Sr = N~^^'^ Y^k^^^^^^k, where the k-sum is over the 
first Brillouin zone of the reciprocal lattice, we obtain from Eq.(|]) for ujn —>■ 

^- J^i:/ws(k,o) , (4) 

where S(k,uj) is the dynamic structure factor of the Cu-spin lattice, 

S{k,uj) = 27rhY,Pa{a\S^\(3)-{(3\S.MS{hu-Ea + Ep) , (5) 

a/3 

and the form factor is 

/(k) = 2 [1 + cos(M] • (6) 

Note that a transition between the two states of the nuclear spin is described by the ladder 
operators /^ = J^ ± il^, which couple to the transverse components S^ = S^ ^ iS^ of 
the electronic spins. However, for temperatures above the Neel temperature the system is 
rotationally invariant in spin space, which enables us to express the NMR rate in terms of the 
rotationally invariant dynamic structure factor (^). This is a very important point, because 
we would like to calculate the dynamic structure factor in Eq. (|) by means of the spin- wave 
expansion. At the first sight it seems that for two-dimensional Heisenberg magnets at finite 
temperatures this approach cannot be justified, because the naive spin- wave expansion is 
based on the assumption of long-range order, which is rigorously known to be absent at any 
finite temperature 0] . For example, an attempt to calculate the staggered magnetization at 
T > via spin- wave theory leads to a infrared divergent integral, signalling the inconsistency 
of the magnon picture in two dimensions 0. The important point is, however, that spin- 
rotationally invariant quantities are free of infrared divergencies, and can be calculated with 
the help the spin-wave expansion even at finite temperatures. For the classical Heisenberg 



model this has been proven by David [|10|, and it is reasonable to assume that it remains 
true for the quantum model. In this work we shall verify explicitly at two-loop order that in 
the rotationally invariant dynamic structure factor (^) all infrared divergencies that appear 
at intermediate stages of the calculation indeed cancel. 



II. THE SPIN- WAVE EXPANSION FOR THE ^^O NMR RATE 

In this section we shall briefly describe our procedure for calculating the dynamic struc- 



ture factor with the help of the Dyson-Maleev formulation of the spin-wave expansion |]TT 
As compared with the more common Holstein-Primakoff formalism |l^, the Dyson-Maleev 
spin-boson mapping has the advantage that the boson representation of the Hamiltonian 
involves only one- and two-body terms, so that beyond the leading order in the spin-wave 
expansion the number of Feynman diagrams generated with the help of the Dyson-Maleev 
transformation is smaller than in the case of the Holstein-Primakoff formalism. Thus, for 
our two-loop calculation of the dynamic structure factor the Dyson-Maleev transformation 
leads to considerable technical simplifications. Of course, final results for physical quantities 
should be identical in both formalisms. 

As usual, we divide the square lattice into two sublattices, labelled A and B, such that the 
nearest neighbors of any given site belong to the other sublattice. The spherical components 
of the spin operators on the A-sublattice are then represented by canonical boson operators 
Or in the following manner. 



(7) 
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(8) 



where br are again canonical boson operators. Because of the two-sublattice structure, it is 
convenient to Fourier transform on both sublattices separately, 

Or = (iV/2)-i/2^V''X , br = {N/2)-'/'J2'^'''"'K , (9) 

q q 

where the prime indicates that the wave- vector sums are over the N/2 points of the reduced 
Brillouin zone, with the wave-vectors q measured with respect to the antiferromagnetic 
ordering wave- vector 11 = [vr/a, vr/a] (see Fig.|D. Because we would like to express all sums 
consistently in terms of the wave- vectors of the reduced Brillouin zone, we rewrite Eq.(P as 

^ = ^E' [/(q)^(q, 0) + /st(q)5st(q, 0)] , (10) 

where 5'st(q, cu) = ^(n -|- q, cu) is the staggered structure factor, and /st(q) = /(H -|- q) is 
the corresponding form factor. Below it will become obvious that at low temperatures the 
inverse thermal de Broglie wavelength acts as an ultraviolet cutoff for the q-sums in Eq.(p!0|). 
Note that the thermal de Broglie wavelength (divided by 2n) of an isotropic antiferromagnet 
can be written as 



he 

Ath = -7^ 



(11) 



where c is the spin-wave velocity. Hence, at low temperatures Ath ^ o-i so that the sums in 



Eq.(lO) are dominated by a small circle in the q-plane with radius X^t^ around the origin. 



Obviously the first term in Eq. ( ]I0| ) represents the contribution from fluctuations of the total 
magnetization. Because these are not critical fluctuations, we expect that the corresponding 
contribution can be calculated perturbatively. From Eq.(^ we see that for small q the form 
factor for l/T^ can be replaced by 



/(q) = 4 + 0(g^) 



(12) 



On the other hand, the second term in Eq.(|TO|) represents for |q| ^ a~^ long wavelength 
critical fluctuations of the staggered magnetization. However, the signature of the critical 
fluctuations contained in S'st(q, 0) is to a large extent removed by the form factor, which 
vanishes for small q as 



hM) = {q.af + 0{qt) 



(13) 



In fact, we shall show shortly that the contributions of both terms on the right-hand side 
of Eq.(|TDp have the same order of magnitude, so that the ^^O relaxation rate is determined 
both by non-critical fluctuations of the total magnetization and by short wavelength anti- 
ferromagnetic fluctuations. In contrast, the corresponding form factor for ^^Cu is unity for 
all wave- vectors [^|-|^, so that the rate 1/Tp" is completely dominated by critical antiferro- 
magnetic fluctuations |^. 

The spin-wave expansion for the dynamic spin-spin correlation functions has been de- 



scribed in an impressive paper by Harris, Kumar, Halperin and Hohenberg [jT3[, so that we 
can be rather brief here and refer the reader to Ref. |T^ for more details. After substituting 
the Dyson-Maleev transformation (0J|) for the spin operators into Eq. (|lD and Fourier trans- 
forming the boson operators as in Eq.(|]), the spin Hamiltonian in D dimensions is mapped 
onto the following boson Hamiltonian 



H 



-NDJS^ + i^Si + H, 



r(4) 
'DM 



(14) 



with 



^Si = 2DJSY^ [a\c 



bifid + 7q(aq^q 



a\b\] 



(15) 



H, 
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N 



1...4 



(16) 



where 1 for qi (and similarly for the other labels), and we have introduced the standard 
notation 7q = -D~^ Si^i cos(q ■ aj). The symbol 5^,(q) indicates momentum conservation 



up to a vector of the reciprocal lattice associated with the sublattices [Tl]. The quadratic 



part -f^£)M of th^ Dyson-Maleev Hamiltonian is now diagonalized by means of a Bogoliubov 
transformation, 





(17) 



with 
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where e^ 



7q 



2 Then we obtain 



H, 



(2) 
DM 



-2DJS 



E'^. 



Q^qClq ~r PqPq 



(19) 



Here iJq = 2D JS'cq is the free magnon dispersion. Taking into account that the /^-operators 
in Eq.(|T9|) are anti-normal ordered, we obtain the following l/S-correction to the ground- 
state energy, 6Eq = —2DJSC, with C = 1 — {2/N)J2q^q- The numerical value of C has 



first been calculated by Anderson W&. In two dimensions the result is C = 0.158. In our 



two-loop calculation presented in Sec.|I^ we shall simply ignore similar zero-temperature 
renormalizations, which involve higher powers of 1/S. Quantum corrections of this type are 
implicitly taken into account by identifying the spin-wave velocity c and the spin stiffness ps 
in our final result (see Eq.(^) below) with the experimentally measured values. Substituting 
the Bogoliubov transformation ( p!7| ) into the quartic part -^dm ^^ obtain totally 10 different 
terms describing scattering of the magnons Oq and /3q in various combinations. The vertices 
for these scattering processes have been derived in Ref. [|T^, and will not be reproduced 
here. For our two-loop calculation we shall use a particular symmetric parameterization of 
the Dyson-Maleev vertices given in Ref. |T^ . It turns out that the complete interaction part 
of the Dyson-Maleev Hamiltonian can be parameterized in terms of five different vertices, 
which can be written as 



l^(^)(q,k,kO = ^'^+^^'^'-V (^)(q,k,kO , J 
MkMk' 



1 5 



(20) 



where the rescaled vertices V^^"^ are non-singular. For our purpose it is sufficient to know 



the leading behavior of these vertices for small wave- vectors |]T3|,|T6|,|T7| 



\/(^-)(q,k,k') 
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(21) 



where k = k/|k| and k' = k'/|k'| are unit vectors. 

The boson representations of the spin-spin correlation functions are obtained in a sim- 
ilar manner. To calculate the dynamic structure factor at finite temperatures, we use the 
fiuctuation dissipation theorem to express the dynamic structure factor in terms of the 
corresponding imaginary frequency spin Green's function. 



2h 

5(q, ^) = [1 + n{hu;/T)] — Im [^(q, luJn)],^^^^^,,^ , (22) 

where n{x) = [e^ — 1]~^ is the Bose-Einstein occupation factor, and 

^(q,za;„) = ^T f'^^ dTe''^^-'"'^"-^^Tr{e-^/^Sr{T) ■ So(0)} . (23) 

J. Jo /j ^ ' 

Here Sr(T) = e^'^Sr(0)e~^'^, and ujn = 2nnT are bosonic Matsubara frequencies. The 
relation between the staggered structure factor S'st(q, uj) and the corresponding spin Green's 
function Qst{q,iUn) is identical with Eq.(|22|). The definition of the staggered spin Green's 
function can be obtained from Eq.(p3D by inserting an additional factor of e*^"" in the sum. 
The spin Green's functions Q{c[,iuJn) and Qst{q, iun) can be calculated via a conventional 
time-ordered perturbation expansion within the Matsubara formalism for bosonic many- 
body systems. 

III. THE LEADING BEHAVIOR OF 1/Tp 
According to Eqs.(p^p^) the ^^O NMR-rate is at low temperatures given by 

(24) 



1 Al 



where 



^F(T) + 1f,,(T) 



nT) = ^Y.'S {ci,0) , (25) 



N\ 



^st(T) = |E'^^st(q,0) . (26) 

From Eqs.(|^) and (^ it is clear that within the Dyson-Maleev transformation one obtains 
in general two distinct contributions to the dynamic structure factor, 

S{q,uj) = Si{(i,uj) + S2iq,uj) , (27) 

(and analogously for S'st(q, uj)), where the first term S'i(q, u) is due to the one-magnon part 
of the spin Green's function Q{q,iu!n), while S'2(q, cu) is due to the two-magnon part (i.e. it 
involves a product of two magnon creation and two annihilation operators). Accordingly, we 
shall refer to 5*1 (q, u) as the one-magnon part, and to 5*2 (q, w) as the two-magnon part of the 
dynamic structure factor. Note that the averages are defined with respect to the interacting 
magnon Hamiltonian, so that both terms involve also higher order virtual magnon excita- 
tions. However, at the one-loop order it is sufficient to evaluate the expectation values with 
the free magnon Hamiltonian -ff^M- ■'■^ ^^^^ approximation the one-magnon term 5*1 (q, a;) 
as well as the transverse part 5*2, ±(q, uj) of the two-magnon contribution vanish in the limit 
cij — i> 0, so that these terms do not contribute to Eq. (pSj) . After a simple calculation we thus 
obtain to leading order 



F(T)^F«(T) = / dqq dkkn{\,^k)[l + n{\,^q)]6{hc{k - q)) , (28) 

TT JO JO 

where the superscript (1) indicates the one-loop approximation. We have written E]f^/T = 
Ath|k|, so that it is obvious that the Bose-Einstein factors act as cutoff functions which 
ehminate the contribution from modes with wavelengths short compared with the thermal 
de Broglie wavelength Ath- Using the fact that by symmetry we may replace g^ — > q^/2 
in Eg. (|26|) , we obtain for the corresponding staggered function FstiT) exactly the same 
expression as in Eq.fpSl). Scaling out the temperature dependence by defining x = qXth, we 
obtain to leading order 

(29) 

O C \ Ath/ <3 C \ ItC / 

Substituting this into Eq. 



3 c VAth/ 
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[he) 


Eq.(24), we finally obtain 18 
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27r 
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2.09 . 



(30) 

We would like to emphasize two points: First of all, from Eq.(p8|) it is obvious that the 
T^-behavior of 1/T° is a simple consequence of the two-dimensional phase space and the 
energy conservation of antiferromagnetic magnons with linear energy dispersion. The second 
important point is that the numerical value of the prefactor in Eq. (|30|) is determined by 
non-critical fluctuations of the total magnetization and by fluctuations of the staggered 
magnetization with typical wavelengths of the order of 1/Ath- From Eqs.(^) and (|28|) it 
is clear that to leading order the relative weight of the total magnetization fluctuations is 
exactly twice as large as the weight of the staggered fluctuations. 

The experimental confirmation of Eq.(pOD would give direct evidence of the existence of 
short wavelength magnons in two-dimensional Heisenberg antiferromagnets at low tempera- 
tures. In the recent experiment by Thurber et al. [|l| the ^^0-NMR rate in Sr2Cu02Cl2 was 
measured with sufficient accuracy to confirm the T^-behavior above the Neel temperature, 
and to obtain an estimate for the numerical value for the prefactor Ci. The experimental 
result is HI Ci ~ 2.7 ± 0.9, where the experimental uncertainty of 30% is mainly due to the 
ambiguities in the measurements of the hyperfine-coupling Aq and the exchange coupling J 
I9| . We conclude that within the experimental accuracy the lowest order spin- wave result 



agrees with the experiment [|l|. This is a very important result, because it demonstrates 
that, to a very good approximation, short wavelength spin- waves in two-dimensional Heisen- 
berg magnets can be treated as free particles. We would like to emphasize that there is no 
long-range order in the system, so that well-defined magnons in two-dimensional Heisen- 
berg magnets at finite temperatures refiect fundamentally different physics than magnons 
in three-dimensional magnets below the Neel temperature. In fact, in two dimensions and 
at T > magnons cannot propagate over length scales larger than the correlation length 
2^ ^ oc exppvrps/T], where ps is the spin stiffness. In other words, at short distances 



the interactions between the magnons become weak, while at long distances the magnons 
completely disappear from the physical spectrum. 



IV. THE TWO-LOOP CORRECTION 



As noticed in Ref. 0, the two- loop correction to Eq.(pUD yields a contribution propor- 
tional to T^, which for an ideal two-dimensional system is negligible at sufficiently low 
temperatures. Although in the appendix of Ref. [^ it was shown how the two-loop cor- 
rection can be obtained in principle, the numerical value of the coefficient of the T'^-term 
was not calculated, because a few years ago sufficiently accurate experimental data to test 
the theoretical prediction were not available. Thurber et al. |I|] have reported for the ffist 
time high-quality data of l/T^ in an almost ideal two-dimensional quantum Heisenberg 
antiferromagnet, which are accurate enough to allow for a detailed comparison with theory. 
Motivated by this experiment, we shall now calculate the numerical value of the two- loop 
correction. 

In quantum antiferromagnets the evaluation higher-order corrections to the leading terms 
in the spin wave-expansion is a rather difficult task, because after the Bogoliubov transforma- 
tion (p^ ) the two-body part ^dm of the spin-wave Hamiltonian (see Eq.(|T6|)) and the boson 
representation of the spin operators involve many different terms. Although the one-loop 
calculation can be performed in a straight-forward manner without resorting to elaborate 
many-body techniques, the classification of the large number of terms contributing to the 
two-loop correction is greatly facilitated with the help of a graphical representation in terms 
of Feynman diagrams |1^ . 



A. The two-magnon part of the total spin structure factor 



We begin with the calculation of the two-loop correction to the two-magnon part 5*2 (q, 0) 
of the total spin structure factor. This is obtained by contracting the interaction part 
-^DM with the two-magnon part of the spin-spin correlation function. Diagrammatically the 
longitudinal and the transverse components of the dynamic structure factor are calculated 
separately. 



5f(q,0) 



4S(q,o) 



4?(q,o) 



(31) 



The superscript (2) denotes the two-loop approximation, while the subscripts 2zz and 2_L 
denote the longitudinal and transverse part of the two-magnon term in the dynamic structure 
factor. Denoting by F^^^ (T) = F\_ (T) -|- FJ^^ (T) the corresponding two- loop correction to 
the dimensionless function F{T) defined in Eq.(^), we have at two- loop order 



F{T) ^ F^^\T) + F^^\T) 



(32) 



where the one-loop contribution F^^\T) is given in Eq. (pQ]) . Using the fact that in two 
dimensions and for large S the spin-wave velocity and the spin-stiffness are related via 
he = 2\/2psa/S, we find after some tedious algebra [|l7l that the transverse part can be 
written as 



Ff^(T) 




Ki 



dx / dyx n{x)n{ii)[l + n{x + y)] 
Jo 



(33) 



and the longitudinal part is 

An /Tn\^ / T \ r /-oo /-oo 

^i'H^^) = -- (^) ( ^ ) [-^2 - i^3 - y^ dxj^ dyxMxMy)[l + nix)] 



(34) 



Here Ki, K2 and K3 are numerical constants of the order of unity, which are given in 
the Appendix. As already mentioned in Sec.|l|, zero temperature l/S'-renormalizations are 
implicitly taken into account in our final result ( [5D| ) by identifying c and ps with the experi- 
mentally measured parameters. The last integrals in the square braces of Eqs.(^) and (|34D 
are logarithmically divergent. This is not surprising, because F]_ (T) and F^^{T) are not 
rotationally invariant, and are therefore sensitive to the fact that our two-dimensional spin 
system does not have long-range order at any T > 0. As discussed at the end of Sec.|, the 
logarithmic divergence should cancel in the rotationally invariant function F\_ {T) + F^f>(T). 



Combining Eqs.(^) and (Q), it is easy to see that this is indeed the case, and we finally 
obtain 

where the finite numerical constant K4 is given in Eq. ( |A6| ) . Comparing Eq.(|35[) with Eq. (p^) , 
we conclude that the two-loop correction involves an additional factor of T/{27Tps), so that 
it is indeed negligible at sufficiently low temperatures. However, as discussed below, for a 
quantitative comparison with the experiment III] this correction can be substantial. 



B. The two magnon part of the staggered structure factor 

According to Eq . (plf) , the rate 1/T^ receives also a contribution from the weighted 
average Fst(T) of the staggered structure factor 5'st(q, 0) defined in Eq.(^). Let us first 
calculate the two- magnon part 5*31,2 (q,i^)- Similar to Eq.(^), we obtain at two- loop order 

>5s?,kq, 0) = 5gL,(q, 0) + ^^^^^(q, 0) . (36) 

The corresponding contributions ivt2±(^) ^^^ -Fst,2zz(^) ^^ the function Fst(T), are 

fSIAT) = F'i'\T) , (37) 

where the numerical constant K^ is given in Eq.( |A7| ), and the divergent integral J is given 
in Eq . (|A9|) . Combining Eqs. (|37D and (^8]), we obtain 

fSim . FSkm . m.m ^ \t (^ [^^ |a, - k. - u . (39, 

where the finite numerical constant K^ is given in Eq. ([All|) . 
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C. The one-magnon part of the staggered structure factor 



There exists one more contribution to the function F^ (T), which is due fact that at 
two- loop order the magnon self-energy acquires a finite imaginary part. This generates a 
finite one-magnon contribution 5'st,i(q, 0) to the staggered structure factor, which in turn 
gives rise to an additional two- loop contribution Fg\. {(T) to Fst(T). This contribution should 
be added to Eq.(p9D in order to obtain the total two- loop correction. Hence, at two-loop 
order we have (compare with Eq.(p^)) 



^(1) 



(2), 



.(2), 



FAT) ^ F^:\T) + F'J!,iT) + F,\^{(T) 



(40) 



By symmetry, no such one-magnon contribution arises in the corresponding total spin struc- 
ture factor. To see this, let us follow Harris et al. [|I^] and denote by ^^^(q, u; -|- iO"*") the 
usual retarded self-energy for a-magnons, and by ^^.^(q, cj -|- iO~^) the off diagonal magnon 
self-energy. The magnon damping can then be written as hT^u{q, u) = — 2ImS^j,(q, cu+iO"''), 
where n^v = a,p. The one-magnon contribution to the total and staggered dynamic struc- 



ture factor at vanishing frequency is |ll3|Jl7 
Si{q,0) = 2hS 



- hm 



q 



where we have defined 

f^^(q,tu) 



Xci 



E^ 
q 



r 



flhl 



lim 



.q>^j 



hVaaiq., UJ) + hTap{q, Uj) 



(41) 
(42) 



'hw/T 



fi,^ = a, P 



(43) 



Noting that |T3[ Taai^i, ^) = ^afsi^., ^^); "we see that for u; ^ the one-magnon contribution 
to the total spin structure factor vanishes by symmetry, 5*1 (q,0) = 0. On the other hand, 
the one-magnon contribution to the staggered structure factor has a finite contribution at 
vanishing frequency, which can be written as 



S,t,i{(i,0) = AhS 



q 



lim hTr 



iq,^j 



(44) 



Calculating the (rescaled) damping rate Taa{^.,0) at two-loop order |T3|JT7|| , we obtain after 
some tedious algebra 



^i(-)-^7(S)'( 



T 



2nps 



with the numerical constant K-/ given in Eq. (|A12|) . 



K-r 



(45) 
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D. The final result 

Combining the results of the previous three subsections with the one-loop result given 
in Eq.(^), we conclude that at low temperatures the two- loop corrections to the functions 

F{T) and Fst{T) can be written as 

F(^\T) = K ( ^] F^'\T) , (46) 

\2TrpsJ 

fS\t) = K,J-^]f,^1\t) , (47) 



with 



K=-^[K^-K2-K,-K^]^ -2.19 , (48) 



^2 

K,, = -^[K^-K,-K, + Kj]^-1.27 . (49) 

From Eq.(^) we thus conclude that at low-temperatures the ^^O NMR rate is given by 

(50) 

If OIL c \iLc/ I \z,nps/ 

with 

C2='^K+^K,,^-im . (51) 

This is our main result. Note that the energy scale of the two- loop correction is set by the 
spin-stiffness p^, so that a measurement of the sub-leading correction can in principle be 
used to obtain the spin stiffness of the material. Because the coefficient C2 is negative, the 
two-loop correction leads to a reduction of the NMR-rate. In other words: Non-interacting 
spin-wave theory over-estimates the magnitude of the NMR-rate of oxygen. Although this 
reduction becomes negligible at sufficiently low temperatures, it is substantial in the exper- 



iment of Thurber et al. |]l|. A simple estimate |19] for the spin stiffness in this experiment 
yields 27rps ^ 1500i^. Taking into account that T k, SOO/T in the low-temperature regime 
above the Neel temperature, we conclude that the correction in Eq.(|50D leads to a reduction 
of the one- loop result by almost 40%. Because the measured rate was approximately 30% 
larger than the theoretical one- loop result, we conclude that the agreement between theory 
and experiment becomes worse if the two-loop correction is included. 

V. SUMMARY AND CONCLUSIONS 

Motivated by recent measurements of the NMR-relaxation rates in the quasi-two- 
dimensional antiferromagnet Sr2Cu02Cl2 by Thurber et al. Jll], we have calculated the 
^"^O NMR-rate above the Neel temperature with the help of the spin-wave expansion at 
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two-loop order. If magnon-magnon interactions are ignored, a calculation based on the two- 
dimensional isotropic Heisenberg antiferromagnet leads to a quantitative agreement with 
the data within the experimental accuracy. We have pointed out that 1/T^ is dominated by 
magnons with wavelengths of the order of the thermal de Broglie wavelength Ath, so that the 
agreement between lowest order spin-wave theory and experiment shows that short wave- 
length magnons are indeed well-defined elementary excitations in two-dimensional Heisen- 
berg magnets at low temperatures. We would like to emphasize that the existence of finite 
temperature magnons in two dimensions reflects fundamentally different physics than in 
three-dimensional systems below the Neel-temperature. 

We have also calculated the leading correction to the T^-behavior due to magnon-magnon 
interactions. This correction involves an additional power of T/{2ttPs), i.e. the energy scale 
for the correction is set by the spin stiffness ps- The numerical value of the prefactor of the 
T^-term is calculated here for the first time. Our main result is given in Eqs. (|50D and (p|), 
and implies that the quantitative agreement between spin-wave theory and the experiment 
becomes less convincing if magnon-magnon interactions are taken into account. It should be 
kept in mind, however, that we have modeled the system by an isotropic two-dimensional 
Heisenberg antiferromagnet. In the low temperature regime close to the Neel temperature 
the various anisotropics and the finite interchain coupling in the experimental system will 
certainly become important, so that Eq.(^) cannot be trusted for temperatures too close to 
the Neel temperature. This might partially explain the slight discrepancy between theory 
and experiment. Another potentially important contribution to the NMR rate, which is not 
contained in our perturbative spin- wave calculation, is due to spin-diffusion 0. However, 
this contribution should be suppressed by anisotropics. Very recently, a measurement of 
the frequency-dependence of l/T^ in Sr2Cu02Cl2 showed that the contribution from spin 
diffusion is indeed negligible |T^. 
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APPENDIX: 

In this appendix we give the expressions for the numerical constants Ki, i = 1, . . . , 7 
introduced in Sec.|rv|. The numerical constant Ki in Eq . (^31) is given by 



1 /"OO /"OO 

-^1 = 0/ ^^ dyxyn{x)n{y)[l + n{x + y)] ^ 1.80 . (Al) 

2 JO Jo 



The constants K2 and K3, that appear in Eq.(p^) can be written as the following four- 
dimensional integrals. 
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K2 = ^ I d^^f dPy xyn{x)n{y) [1 + n{x)] , (A2) 



K, = —ld'^d'y-n{x)[l+n{x)y-^ ^^ . (A3) 



Here q is an arbitrary fixed unit vector. The integral ([A2D can be reduced to a two 



dimensional one with the help of circular coordinates in the x- and y-planes, which is then 
easily calculated numerically, 



oo 



K2= dx dy xyn{x)n{y)[l + n{x)] ^ 2m . (A4) 

JO Jo 

The integral K^ can be calculated analytically by shifting y — »• y — r/2, with r = x — xq, and 
then introducing circular coordinates (x, y?) in the x-plane, and elliptic coordinates (m, a) 
in the y-plane, i.e. y ■ f = | cosh m cos a, and y • f_|_ = |sinhMsina. Here r = |r| = 



xJ2{l — q-x), and r_L is a unit vector perpendicular to f. With these coordinates the 
integrations in K^ factorize and can then be done exactly, with the result 

^3 = ^C(3) ^ 1.80 . (A5) 

The constant K^ in Eg . (|35|) is given by 

/■oo /'OO 

K4 = dx dyx^n{x)n{y)[n{x) ~ n{x + y)] ~ 1.52 . (A6) 



Note that the first term in the square brace of Eq. ( |A6|) is due to longitudinal fluctuations, 
while the second term is generated by transverse fluctuations. Each of these terms sepa- 
rately would lead to a logarithmically divergent integral, but the divergence cancels in the 
rotationally invariant function F^'^^T), in agreement with Ref. fl^. This is also a non-trivial 
check that in our two-loop calculation we did not forget any diagram. 
The numerical constant Kr, in Eq.(^) is 

A-,^--L/d^x//y^nM|l + nM| '.^-*-'''*-'"q-"-^> , . (A7) 

47r2 y J •" y ^ '^ ^ ^^(|a;q-x-y| -F7/)|a;q-x-y| ^ ^ 

This is a constant of the order of unity, which can be calculated analytically with the help of 
circular coordinates in the x-plane and (after shifting y -^ y — |(x — xq)) elliptic coordinates 
in the y-plane. We obtain 

K5 = 4C(3) ^ 4.80 . (A8) 

The integral J in Eq.(^) is 



1 f£^ fd^y^nix)[l+nix)Hy)—^^—^^ ^ 
2tt^ J J y |xq — X — yp — ?/"^ 



J=— I £^l d^y-nix)[l+nix)]niy)j-^;^---^^—^ . (A9) 



Introducing circular coordinates in the x- and y-planes, the two angular integrations can be 
performed analytically, with the result 
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J 



dx dy{x'^ + y'^)n{x)[l + n{x)]n{y) . 



(AlO) 



Note that for small y the Bose-Einstein factor n(y) diverges as 1/y, leading to a logarithmic 
divergence of the integral J. However, in the rotationally invariant correlation function this 
divergence is again cancelled, so that the constant Kq in Eq.(P^ is finite. 



K. 



oo rx 

dx / dy{x'^ + y'^)n{x)n{y)[n{x) — n{x + y)] ~ 0.89 . 

JQ 



(All) 



Finally, the numerical constant K^ in Eq.(HS|) is given by 



1 /"OO /"OO 

-^7=0/ ^^ I dyx'^n{x)n{y)\i + n{x + y)] 
2 Jo Jo 

^ 1.80 . 



1 + 



X + y 



X 



^x"^ + Ay{x + y) 



(A12) 
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FIGURES 




FIG. 1. Cu-0 planes and hyperfine couplings. The black dots mark the positions of the 
oxygen atoms, and the empty circles are the copper atoms with the electronic spins. The shaded 
areas symbolize the hyperfine coupling between a given oxygen site and the two neighboring copper 
spins. 




FIG. 2. First Brillouin zone of a square lattice and reduced Brillouin zone of the antiferro- 
magnet (shaded area). In the reduced zone scheme all wave- vectors q are measured with respect 
to the antiferromagnetic ordering wave-vector 11 = (7r/a,7r/a). 



